Mean-field instability of trapped dilute boson-fermion mixtures 
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The influence of boson-boson and boson-fermion interactions on the stability of a binary mix- 
ture of bosonic and fermionic atoms is investigated. The density profiles of the trapped mixture 
are obtained from direct numerical solution of a modified Gross-Pitaevskii equation that is self- 
consistently coupled to the mean-field generated by the interaction with the fermionic species. The 
fermions which in turn feel the mean-field created by the bosons are treated in Thomas-Fermi ap- 
proximation. We study the effects of different combinations of signs of the boson-boson and the 
boson-fermion scattering lengths and determine explicit expressions for critical particle numbers as 
function of these scattering lengths. 

PACS numbers: 03.75.-b, 03.75.Fi, 32.80.Pj 



Recent experimental successes in the trapping and 
cooling of mixtures of bosonic and fermionic atoms 
|[]. H constitute a new branch in the field of trapped ul- 
tracold gases. Similar to the purely bosonic gases boson- 
fermion mixtures offer unique possibilities to study fun- 
damental quantum phenomena. Moreover they appear 
as a promising candidate to realize a BCS transition to 
a superfluid phase of the fermionic component One 
of the most appealing features of these systems is that 
the strength of the interaction between the atoms can be 
tuned in a wide range by utilizing a Feshbach resonance 

For the sympathetic cooling of a Fermi gas in binary 
boson-fermion mixtures the collapse caused by attractive 
interactions is responsible for a severe limitation of the 
lowest achievable temperature |l], || . The occurrence of a 
mean-field instability when the density or particle num- 
ber exceeds a critical value was already studied experi- 
mentally for purely bosonic systems || . In this communi- 
cation we investigate the interplay between boson-boson 
and boson-fermion interactions and the implications for 
the stability of the mixture. 

In order to describe the properties of the binary boson- 
fermion mixture at zero temperature we first construct 
the energy functional in mean-field approximation. The 
many-body state |\&) of the mixed boson-fermion system 
is a direct product of a symmetric -/Ve-body state |\&b) 
for the bosonic species and an antisymmetric -/Vp-body 
state I^f) for the fermions. The Hamiltonian of the 
interacting mixture reads 
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The first line contains the kinetic energy operators and 
the external trapping potentials U B (x) and Up(x) for the 
bosonic and the fermionic species with masses tub and 
tof, respectively. 



Due to the large average distance the atom-atom in- 
teraction can in general be described by an effective con- 
tact interaction for all partial waves |Q ^]. In a binary 
boson-fermion mixture one has to distinguish three in- 
teraction types: boson-boson interactions, boson-fermion 
interactions, and fermion-fermion interactions. The s- 
wave interaction between two bosons is described by the 
third term of Eq. (|l]), where ob is the s-wave scattering 
length. Since we consider a pure Bose-Einstein conden- 
sate at zero temperature only the s-wave term is needed. 
Higher even partial waves are negligible. For the inter- 
action between a bosonic and a fermionic atom s- and 
p-wave terms contribute. The operator of the s-wave 
boson-fermion contact interaction forms the last term of 
the Hamiltonian (|j), where obf is the corresponding s- 
wave scattering length and tobf = 277Zb^f /(^b + %) 
is twice the reduced mass of the pair. Since the s-wave 
interaction dominates in many cases of interest, e.g. the 
6 Li/ 7 Li mixture, we will neglect the p-wave interaction 
for this discussion. For identical fermions s-wave contact 
interactions are prohibited by the Pauli principle. The 
first nonvanishing contribution is the nonlocal p-wave in- 
teraction, which will also be neglected in the following. 
However, p-wave interactions can have significant influ- 
ence on the structure and stability of the fermionic com- 
ponent as we discussed earlier 0, |[. 

The expectation value of the Hamiltonian ([l]) calcu- 
lated with the many-body state |*) = |* B ) ® |*f), 
defines the energy density of the mixture 



E= (9\H\^f) = J d 3 x £[n B ,n F ](x). 



(2) 



The energy density £ [tib, Tip] is decomposed into a purely 
bosonic part (B), a fermionic part (F) and the interaction 
part between the two species (BF) 

£[riB,n F ] = £ B [n B ] +£f[«f] + £bf [tib , rip] . (3) 

For the calculation of £b["b] we assume that the 
bosons are in a pure Bose-Einstein condensate, i.e. the 
A^B-boson state I^b) is given by a direct product of iden- 
tical single particle states \4>b)- This immediately leads 
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FIG. 1: Radial density pro- 
files of a boson-fermion mixture 
with iV P = iV B = 10 000 for 
different interaction strengths. 
The boson density ran (a?) is 
given by the solid line (left 
scale) and the fermion density 
np(x) by the dashed line (right 
scale) . The upper row shows ex- 
amples with increasing boson- 
fermion attraction and fixed 
a B /e = 0.001. The lower row 
depicts examples with increas- 
ing boson-fermion repulsion and 
a B /e= -0.00004. 



to the standard Gross-Pitaevskii energy density ]Io| 
£b[«b](^) = U B {x)n B {x) + — !— |Vn B /2 (f)| 2 
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where ub(x) — $ B (x) = -/Vb(x|<^>b) 2 is the ground-state 
boson density. 

For the evaluation of the parts of the energy density 
that involve the fermionic species we employ the Thomas- 
Fermi approximation. It was shown that this is an excel- 
lent approximation for the particle numbers considered 
here The fermionic part of the energy density in 

Thomas-Fermi approximation reads 

3 5/3 n A/3 

£ F [n F ](x) = U F (x)n F (x) + , n F ' (x). (5) 

52 L ' A m F 

The interaction between the bosonic and the fermionic 
species yields the contribution 

£BF[nB,n F ](x) — 47raBF n B {x)n F (x). (6) 
msF 

The functional minimization of the total energy (^) 
under the constraint of given numbers of bosons and 
fermions leads to the density profiles of the trapped 
gas. The constraints are implemented by introducing 
the chemical potentials /iB and (i F and minimizing the 
transformed energy 

F = J d 3 x [£{nB,n F }(x) - fiBn B (x) - fJ, F n F (x)] . (7) 

Variation with respect to the fermion density n F (x) im- 
mediately leads to an equation for the density distribu- 
tion 
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From the functional variation of the energy (0) with re- 
spect to the bosonic density we obta in an E uler-Lagrange 
differential equation for $b (x) — y^nB(x) 
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$b(x) = a«b$b(x). 



This is a modified Gross-Pitaevskii equation which in- 
cludes the mean-field contribution generated by the in- 
teraction with the fermionic species. 

The simultaneous solution of the coupled Eqs. (||) and 
(0) gives the density profiles of the two species. We solve 
the nonlinear differential equation (^J) with an efficient 
quantum diffusion algorithm using a fast Fourier trans- 
formation Jljj . The ground state solution is extracted by 
successive application of the time evolution operator for 
a small imaginary time step. We accomplish the solution 
of the coupled problem by a simple iterative procedure: 
© the boson density is initialized with the Gaussian pro- 
file of the noninteracting Bose gas with Nb particles. 
The fermion density is calculated using (JsJ) with /iF ad- 
justed such that the integral over n F (x) gives the desired 
particle number N F . (2) A single imaginary time-step is 
performed using the mean-field according to the boson 
and fermion densities obtained in the previous two steps. 
The resulting hb{x) is normalized to Nb and used as 
initialization for the next iteration cycle. 

With these tools we investigate the instability of the 
boson-fermion mixture against collapse induced by at- 
tractive boson-boson or fermion-boson interactions. In 
order to keep the discussion simple we restrict ourselves 
to spherical symmetric systems with equal numbers of 
bosons and fermions Nb = N F . We assume equal masses 
for the two species m = ttjb = m F = ttzbf and iden- 
tical parabolic trapping potentials U B (x) — U F (x) = 
x 2 /(2m£ 4 ). The oscillator length i — (moj) -1 / 2 serves 
as fundamental length unit for the numerical treatment. 
A more general treatment including asymmetric boson- 
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fermion mixtures in deformed traps will be presented in 
a subsequent paper. 

First we consider the case of repulsive boson-boson and 
attractive boson- fermion interactions (ob > 0, obf < 0). 
Here the attractive interaction between the species in- 
duces a mean-field collapse if the densities or particle 
numbers exceed a critical value. The upper row of Fig. 
|l] shows the density profiles of (meta)stable configura- 
tions with Nb = Np — 10 4 particles for three different 
values of the boson-fermion scattering length obf and 
fixed a^/i = 0.001 (corresponds to ob ~ 20aB O hr for a 
typical trap with I = 1 \xva) . Due to the Pauli principle 
the fermionic density distribution is much more spread 
out and has a significantly lower central density than the 
bosonic distribution with the same particle number (no- 
tice the different scales for ub{x) and tif(x) in Fig. [l]). 
Attractive boson-fermion interactions generate an attrac- 
tive mean-field for bosons proportional to the density of 
the fermions and vice versa. This causes an increase of 
both densities in the overlap region as can be seen in 
panel (b) of Fig. |. With increasing strength of the 
boson-fermion attraction the fermion density grows sub- 
stantially. As shown in Fig. [l](c) the fermion density 
can easily be increased by a factor 3 compared to the 
noninteracting case. 

If the strength of the boson-fermion interaction exceeds 
a critical value then the mixture collapses towards high 
densities. In this case the attractive mean-field is not 
stabilized by the positive kinetic energy contribution or 
the repulsive boson-boson interaction any more, i.e. the 
gas can lower its energy by contracting and increasing 
the density in the central region. This phenomenon was 
investigated in detail for purely bosonic [ [b3| and purely 
fermionic jj], |i| systems. To our knowledge the collapse 
of boson-fermion mixtures was studied only inchoately 
using parameterized density profiles [l4| . 

We can determine the critical boson-fermion scatter- 
ing length obf for which the collapse occurs with a rather 
simple numerical procedure. The quantum diffusion algo- 
rithm used to obtain the solution of the Gross-Pitaevskii 
equation (|^) diverges if the mean-field instability occurs, 
i.e., the change of the boson density in the trap center 
increases for successive imaginary time steps. Thus by 
observing the convergence behavior of the central den- 
sity during the imaginary time evolution we can decide 
whether the mixture is stable or collapses. 

To obtain a simple measure for the stability we proceed 
in two steps: First we determine numerically the critical 
boson-fermion scattering length for a set of particle num- 
bers Np = Nb — 500 . . . 10 6 and boson-boson scattering 
lengths o,b/£ = . . . 0.003. Then we fit a parametrization 
which connects the particle number with the two scatter- 
ing lengths to this data set. This leads to an expression 
for the critical particle number iV cr as function of the 
scattering lengths qb/£ > and obf/^ < 0: 
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FIG. 2: Logarithm of the critical particle number as func- 
tion of the boson-fermion scattering length <ibf 1 1 for differ- 
ent values of the boson-boson scattering length. The solid 
line corresponds to a^jl = 0. Thick lines show the behavior 
for repulsive boson-boson interactions: a-^ji — 0.0003 (thick 
dash-dotted), 0.001 (thick dashed) , 0.003 (thick dotted) . Thin 
lines correspond to boson-boson attraction: as/i = —0.0001 
(thin dash-dotted), -0.0003 (thin dashed), -0.001 (thin dot- 
ted). 



Any mixture with Nb = Np > N CI is unstable against 
mean-field induced collapse. 

Figure || depicts the critical particle number as func- 
tion of the boson-fermion scattering length obf for dif- 
ferent values of the boson-boson scattering length as- 
The thick lines show the critical particle number obtained 
from (^|) for repulsive boson-boson interactions. Notice 
that the values of cib/£ used in the plot are typically 
one order of magnitude smaller than the range shown for 
clbf/(', this emphasizes the strong influence of the boson- 
boson interaction. 

We observe that a moderate boson-fermion attraction 
causes a severe limitation of the particle number of the 
(meta)stable mixture. For as/i — and obf/^ = —0.01, 
which corresponds to asF ~ — 200aB O hr for £ — Xfim, 
the mixture is stable only up to TVb = Np w 1000. The 
inclusion of a repulsive boson-boson interaction leads to 
a significant stabilization, i.e. an increase of the critical 
particle number. For the above example, the critical par- 
ticle number grows from iV cr w 1000 to 8500 if a weak 
boson-boson repulsion with ob/^ = 0.001 is included. 

As a second class of systems we consider mixtures 
where both, the boson-fermion and the boson-boson in- 
teraction are attractive (ob < 0, asF < 0). Compared 
to the previous class the attractive boson-boson interac- 
tion enhances the attractive mean-field generated by the 
boson-fermion interaction and abets the instability. The 
critical particle numbers are significantly reduced. 



We obtain a relation between the scattering lengths as 
and obf and the critical particle number in the same way 
as before. The resulting parametrization of the critical 
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particle number for a B < and qbf < reads 



iYcr(aB,<2BF,f) = 



0.575 



|aB/^| + 2.03|a B F/^| 



1.78 ' 



(11) 



For cibf = we recover the well known relation between 
the critical particle number and |as| for a pure Bose gas 
|To| , |l3| . The additional term in the denominator de- 
scribes the modification of iV cr in the presence of an at- 
tractive boson-fermion interaction. In the limiting case 
a B = this relation coincides with the corresponding 
limit of ©. 

The thin curves in Fig. ^ show the behavior of the crit- 
ical particle number for a B /£ < 0. The presence of the at- 
tractive boson-boson interaction reduces the critical par- 
ticle number for a B -p < significantly. With increasing 
strength of the boson-boson attraction the influence of 
the boson-fermion attraction is largely reduced. 

Finally we consider mixtures with attractive boson- 
boson and repulsive boson-fermion interactions (<2b < 
0, cl B f > 0). The 6 Li/ 7 Li mixture used in the experiment 
of Truscott et al. Q belongs to this class of interactions. 
The lower row of Fig. [I] shows the density profiles for 
three different values of a B p > 0. Already for a very 
weak boson-fermion repulsion the two species separate 
spatially [see Fig. (j](e)], the bosons occupy the central 
region of the trap (boson core) and the fermions consti- 
tute a shell around it Jll], |l5). This structure may have 
interesting implications for the mean-field instability of 
the bosons: The fermionic shell compresses the boson 
core, i.e. increases the maximum boson density as it is 
clearly seen in Fig. |l|(f). This could promote the mean- 
field collapse in the presence of attractive boson-boson 
interactions and lower the critical particle number. 

The dependence of iV cr on the scattering lengths, which 
is obtained from the direct numerical solution of the cou- 



pled problem, is for ob < and a B p > very well de- 
scribed by the parametrization 



N CI (a B ,a BF ,£) 



0.575 - 0.230 (obf/£) 

\a B /e\ 



0.333 



(12) 



The thin curves for <i B f/£ > in Fig. g show this 
dependence. Obviously the influence of the repulsive 
boson-fermion interaction on the critical particle number 
is marginal. The critical particle number reduces slightly 
if a B F / 1 is increased. This can be attributed to the com- 
pression of the boson core mentioned before. Although 
the boson-fermion interaction has a strong influence on 
the density profiles, its influence on the critical particle 
number is negligible. 

In summary we have investigated the mean-field insta- 
bility of binary boson-fermion mixtures with equal par- 
ticle numbers. We solved the coupled Gross-Pitaevskii 
equation numerically and obtained the critical particle 
number as function of the boson-boson and the boson- 
fermion scattering lengths. We have shown that the 
boson-boson and the boson-fermion interaction have very 
different effects on the collapse: In the presence of at- 
tractive boson-fermion interactions (a B F < 0) the sys- 
tem is stabilized by weak repulsive and destabilized by 
weak attractive boson-boson interactions. In contrast 
repulsive boson-fermion interactions (obf > 0) cannot 
stabilize the mixture against collapse due to attractive 
boson-boson interactions. With regard to the recent ex- 
periments using mixtures of 6 Li (F — 3/2, mp = 3/2) 
and 7 Li (F = 2, mp = 2) |Q, this implies that the 
collapse due to the attractive boson-boson interaction is 
influenced only marginally by the boson-fermion repul- 
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